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Abstract. Quantum criticality in systems of local moments interacting with itin¬ 
erant electrons has become an important and diverse field of research. Here we re¬ 
view recent results which concern (a) quantum phase transitions in single-impurity 
Kondo and Anderson models and (b) quantum phase transitions in heavy-fermion 
lattice models which involve critical quasiparticles. For (a) the focus will be on im¬ 
purity models with a pseudogapped host density of states and their applications, 
e.g., in graphene and other Dirac materials, while (b) is devoted to strong-coupling 
behavior near antiferromagnetic quantum phase transitions, with potential appli¬ 
cations in a variety of heavy-fermion metals. 


1 Introduction 

Quantum phase transitions (QPT) |1I2I3I4I5IB] continue to be an exciting field of research 
in condensed matter physics, both theoretically and experimentally. While for the simplest 
situations, such as magnetic transitions in Mott insulators devoid of quenched disorder, agree¬ 
ment between theory and experiment has been achieved, more complicated cases are far from 
being fully understood. These include (i) quantum phase transitions with conventional, i.e., 
symmetry-breaking, order parameters in metallic systems [3], (ii) quantum phase transitions 
under the influence of quenched disorder [7], (iii) interaction-driven metal-to-insulator transi¬ 
tions [8], (iv) topological transitions in interacting systems and (v) boundary quantum 

phase transitions m- Each of these cases poses its specific challenges to theory, and the lack 
of a complete theoretical framework prohibits an efficient classification and understanding of 
experimental observations. 

In this paper, we will provide a partial review on strands of work performed in the context of 
Kondo systems [T3] which belong to the cases (i) and (v), namely antiferromagnetic transitions 
in strongly interacting metals and transitions in variants of the single-impurity Kondo model, 
describing isolated magnetic moments immersed in baths of fermionic quasiparticles. The con¬ 
ceptual link is given by the occurrence of “critical” fermionic quasiparticles, where fermionic 
spectral functions display power-law singularities. In addition, we will highlight connections to 
case (iii), namely to Mott transitions and partial (or orbital-selective) Mott transitions where 
critical quasiparticles are expected as well. 

The body of the paper is organized as follows: In Section we will discuss conceptual as¬ 
pects of critical quasiparticles and provide a quick survey on available results in the literature. 
Section is devoted to the Kondo model with a power-law density of states. We summarize 
the theoretical understanding of its quantum phase transitions and discuss some recent appli¬ 
cations, including Kondo impurities in graphene and on the surface of topological insulators. 
Section [previews a recently developed theory for antiferromagnetic bulk quantum phase transi- 
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tions, which assumes strong coupling between fermions and order-parameter fluctuations such 
that power-law singularities develop in the fermionic spectrum on the entire Fermi surface. A 
summary and outlook will be given in Section 


2 Critical quasiparticles 

Phase transitions into symmetry-breaking phases are characterized by the critical behavior of a 
collective order parameter [1]. This order parameter is a bosonic field; for transitions of electrons 
in solids it can usually be represented as a composite of two (or four) fermion operators. In 
contrast, single-particle excitations are not generically critical at the transition point. However, 
recent research has highlighted the possibility - and potential experimental relevance - of critical 
single-particle excitations at certain quantum phase transitions - these form the subject of this 
review. 

In both lattice and continuum systems, power-law behavior of fermionic single-particle prop¬ 
agators may occur as function of frequency or wavevector, or both. As fermionic low-energy 
excitations of stable metals occur on Fermi surfaces (more generally hypersurfaces) in momen¬ 
tum space, critical power-law singularities may similarly occur on surfaces in momentum space: 
Most naturally, the Fermi surface of a metal may develop into a surface of critical quasiparti¬ 
cles upon approaching the quantum critical point This is distinct from the standard bosonic 
order-parameter situation, where singular behavior of correlation functions at criticality is usu¬ 
ally restricted to the vicinity of a single point in energy-momentum space. 

Critical fermionic quasiparticles have been discussed in various contexts. On the one hand, 
they may occur as “active” critical degrees of freedom at Mott metal-insulator transitions: 
This has been proposed within effective theories for Mott and orbital-selective Mott transitions 
|13 | 14 | 15| . For such transitions there is no local order parameter, but the loss of metallicity 
accompanied by a change in the Fermi volume is a defining criterion. Available analytical theo¬ 
ries have employed representations in terms of slave particles and gauge fields; in this language 
the Mott transition is often marked by the condensation of an auxiliary boson coupled to a 
gauge field. Recent numerical findings |I6j of quantum critical scaling of the conductivity near 
the Mott critical endpoint as described by dynamical mean-field theory also hint at fermionic 
criticality. 

On the other hand, fermions of a metal may couple to critical order-parameter degrees 
of freedom in such a way that quasiparticles are rendered critical by this coupling [3]. In 
general, this may happen in a momentum-selective fashion, such that, e.g., only fermions at so- 
called hotspots (momenta connected by the wavevector of the order parameter) become critical. 
More interesting are situations where fermions at all (or almost all) Fermi-surface momenta 
become critical. The latter applies to the particular strong-coupling theory of antiferromagnetic 
quantum phase transitions which will be reviewed in Section It applies as well to Ising- 
nematic quantum phase transitions in 2d m where the ordering wavevector is zero and hence 
all momenta are hotspots. Finally, singular behavior of fermions at all momenta also underlies 
the idea of local quantum criticality in the framework of extended dynamical mean-field theory 

m- 

For quantum phase transitions in fermionic quantum impurity models m. power-law be¬ 
havior of local propagators is not uncommon, and we will review a well-studied case in Sectionj^ 
In some of the models, fermions can indeed be interpreted as the critical degrees of freedom. 
An understanding of these types of impurity criticality (being often more tractable than lattice 
cases) may also be relevant to the critical behavior of lattice models as described by dynamical 
mean-field theory |19j and extensions thereof. 


^ We use the term “critical” in the context of fermionic spectral functions if the quasiparticle weight 
vanishes, i.e., the S peak is replaced by a critical power-law continuum as function of frequency. We 
note that such power-law spectral functions are realized, e.g., in one-dimensional interacting electron 
systems described by Luttinger-liquid theory. 
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3 The pseudogap Kondo effect 

As a well-understood example for critical quasiparticles, we discuss the rich physics of the 
so-called pseudogap Kondo model. This model has appeared first in the context of magnetic 
impurities in unconventional superconductors |20| . While this initial work only uncovered the 
existence of a quantum phase transition between a screened and an unscreened impurity mo¬ 
ment, detailed numerical studies |2ll22l23l24| using Wilson’s numerical renormalization group 
(NRG) technique |25l2(i) determined the complete phase diagram and some of the critical prop¬ 
erties. A complete analytical understanding of the low-energy physics was achieved 
after realizing that the Anderson model provides a framework which allows to utilize controlled 
epsilon-expansion techniques. Later, the pseudogap Kondo model was applied to describe im¬ 
purities in graphene and other Dirac materials |29| . 

In the following we review the basic aspects of the pseudogap Kondo problem, with an eye 
on critical quasiparticles, and discuss recent applications. We will restrict our attention to the 
case of a spin S' = 1/2 coupled to a single screening channel; the two-channel version has been 
discussed in detail in Refs. HIED]. 


3.1 The pseudogap Kondo and Anderson models 


The standard Kondo Hamiltonian [12] reads 

■H = ^ ekcLck^ + yo cLck'<r + JoS- Sq, (1) 

k,(T k.k',(7 

where S is the impurity spin S = 1/2, Cko- are conduction-electron operators, and Sq = 
I Skk' '^k(T'^'^<T"'k'cr' their spin density at the impurity site, with r the vector of Pauli matri¬ 
ces. The Kondo coupling Jq and the potential-scattering strength Vq characterize the impurity. 
In the following, we denote the density of states (DOS) of the conduction band at the impurity 
site by p(a;). 

The Kondo model may be derived from the more general Anderson model, 

n = E ekcLck<T + + UUdtndi + ^ (vkcl^d^ + h.c^ (2) 

k,(T C7 k,(T 


which describes an impurity level d, with Uda = d]yd„, hybridized with a conduction band. In 
this model, the d level prefers single occupancy if < 0 and Cd + U >0. For large |ed| and t/, 
charge fluctuations are frozen out, such that an effective spin 1/2 degree of freedom remains 
whose dynamics can be mapped to a Kondo model ffl . Second-order perturbation theory yields 

m- 


Jo = 2v^\^ 
\^d\ 


\u- 


■edl 


Fn = ^ ( ^ - 


\U -\- Cd\ 


(3) 


where Uk = u has been assumed. 

For a metallic host, the DOS p(uj) is finite at the Fermi level. Then, for antiferromagnetic 
Jo > 0 the impurity spin is screened below the so-called Kondo temperature Tk- For a flat 
conduction-band DOS, p{uj) = Po, one finds [ID] : 


Tk = V^i^e-F(APo) . (4) 

In the following, we will instead concentrate on the case of Kondo impurities coupled to 
fermions with a pseudogap DOS, of the low-energy form p{(jj) cx For r > 0 this DOS 
corresponds to a semimetal; it can be thought of arising from linearly dispersing Dirac electrons 
in (1 -|- r) space dimensions, with r = I corresponding to graphene. As a consequence of the 
vanishing DOS at the Fermi level the tendency toward Kondo screening is reduced, such that no 
screening occurs at small Kondo coupling Jq, and a quantum phase transition between phases 
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without and with screening occurs upon increasing Jq [ni2o] . As the critical behavior discussed 
below is independent on high-energy details of the DOS, we parameterize p{uj) as 


p{u:) 


I + r 
2D^ 


\ujrei\oj\-D) 


(5) 


with D being a bandwidth. 


3.2 Phase diagram 


Despite its simplicity, the pseudogap Kondo model has an extraordinarily rich phase diagram, 
first determined by Gonzalez-Buxton and Ingersent [53] using NRG. The physics depends not 
only on Jq and the exponent r of the low-energy DOS, but also on the presence or absence of 
particle-hole (p-h) symmetry. Importantly, p-h symmetry requires both p(pj) = p{—uj) in the 
host and Vq = 0 in the Kondo model Q [or 17 = —2ed in the Anderson model (2|]. 

We start by summarizing the phase diagram of the pseudogap Kondo model. We restrict 
our attention to the case of antiferromagnetic Kondo coupling, Jq > 0, and use the acronyms 
of Ref. |53| for both phases and fixed points. 

The phase diagram, qualitatively sketched in Fig. [l] shows that Kondo screening is only 
realized for large Jq; in addition p-h asymmetry is favorable for screening. There are three stable 
phases: the local-moment phase (LM), the p-h symmetric strong-coupling phase (SSG), and the 
p-h asymmetric strong-coupling phase (ASC). Here, LM corresponds to an asymptotically free 
(unscreened) spin, with a residual entropy 5'imp = ln2. SSG represents the generalization of the 
metallic Kondo-screened phase to finite r. Due to the vanishing host DOS, SSG is characterized 
by a residual entropy S’imp = 2rln2, implying partial screening. Moreover, SSG only exists 
for r < 1/2. Finally, ASC is a fully screened phase with vanishing residual entropy. In stark 
contrast to the metallic case r = 0 where p-h asymmetry is marginally irrelevant, we see that 
p-h asymmetry is relevant for r > 0, i.e., SSG is destabilized by any amount of p-h asymmetry. 

As deduced from the numerical solution of the pseudogap Kondo model |53|, the topology 
of the phase diagram changes qualitatively as the bath exponent r is varied, see Fig. [l] 


a) 0 < r < Tniax = 1/2. For particle-hole symmetry, a critical coupling Jc separates LM from 
SSG: For initial values J < the system is in the LM phase, whereas it is in the SSG phase 
for J > Jc- For finite particle-hole asymmetry, i.e. Vq ^ 0, there is a boundary between LM 
and ASC. 


b) r > Tmax- SSG disappears such that there is no Kondo screening at particle-hole symmetry, 
irrespective of the strength of the Kondo coupling Jq. In contrast, screening is still possible for 
finite asymmetry, where a transition between LM and ASC continues to exist. 


c) — 1 < r < 0. This regime can possibly be realized in the case of reconstructed vacancies 
in graphene [3TI35I but was analyzed more generally in Ref. [33]. SSG is stable, and a phase 
transition separates SSG from a newly emerging phase ALM, corresponding to Jq = 0 and large 
particle-hole asymmetry. In the following we will, however, not discuss r < 0 in any detail. 


3.3 Quantum phase transitions 

The pseudogap Kondo model features two qualitatively different transitions: The transition 
between LM and SSG is controlled by a p-h symmetric critical (SCR) fixed point, while the 
transition between LM and ASC is controlled by a p-h asymmetric critical (ACR) fixed point 
(for r > r*, see below). The existence of SCR is tied to that of the SSG phase; SCR is thus 
present for 0 < r < 1/2. 
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a) 0 < r < = 1/2 



b) r > r, 

Voi 



ASC 

LM 



0 - J, 


Fig. 1. Schematic phase diagram for the pseudogap Kondo model in the plane spanned by the Kondo 
coupling Jo and the potential scattering Vo, the latter measuring particle-hole asymmetry. The phases 
are denoted as local-moment phase (LM), symmetric strong-coupling phase (SSC), and asymmetric 
strong-coupling phase (ASC). For 0 < r < Tmax = 1/2 screening can occur at particle-hole symmetry 
(left), whereas screening is only possible in the presence of particle-hole asymmetry for r > rmax (right). 
In addition to the change at Tmax, the structure and nature of the critical fixed points also changes at 
r* = 0.375 ± 0.002 and r = 1, for details see text. 


Interestingly, also ACR exists in a restricted r range only, namely for r > r* k, 0.375. As 
r approaches r* from above, ACR merges with SCR, such that for 0 < r < r* the transition 
between LM and ASC is controlled by SCR which is then a multicritical point. Further, the 
critical exponents of ACR become trivial for r > 1, showing that r = 1 plays the role of an 
upper critical dimension |27j . 

The pseudogap Kondo and Anderson models share identical hxed points and quantum phase 
transitions [53]. This observation can be rationalized within the effective field theories described 
in what follows. 


3.4 Critical field theories 

The complicated topology of the RG flow |23| suggests that different field theories are required 
to describe the critical properties near the SCR and ACR hxed points. Such held theories 
have been worked out in detail in Refs. [57153] and provide an essentially complete analytical 
understanding of the pseudogap Anderson and Kondo models. Interestingly, none of these 
held theories is of conventional (i.e. bosonic) Landau-Ginzburg-Wilson type; instead all are of 
genuinely fermionic character and are formulated in the degrees of freedom of either the Kondo 
or the Anderson model. 

In the following we shall summarize the three relevant critical theories. When specifying 
how equations from perturbative RG, we will assume a symmetric pseudogap density of states 
as in Eq. (§. The effect of a high-energy particle-hole asymmetry in the DOS can be absorbed 
in the impurity part of the Hamiltonian, e.g., the potential scattering term of the Kondo model. 
This can be rationalized within RG, where integrating out the particle-hole asymmetric piece 
of the bath at high energies yields an effective model with particle-hole symmetric bath at 
low energies and a renormalized impurity Hamiltonian, where in particular the particle-hole 
asymmetry is accumulated. 


3.4.1 SCR: Kondo model 

For small r, an efficient description of the physics at SGR is obtained via the Kondo model itself, 
Eq. ([^. A perturbative expansion can be performed in Jq and Vq around the LM fixed point 
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where Jq = Vq = 0 [20124128] - this can be understood as a generalization of Anderson poor 
man’s scaling [33] to the pseudogap case. As is standard practice, we introduce dimensionless 
couplings j and v, for details see Ref. jSSj. Power counting reveals that both couplings are 
marginal for r = 0 and irrelevant for r > 0, dim[j] = dim[u] = —r. The one-loop flow equations 
read 


dj 

d\a.D 


= rj - f 


and 


dv 

dlnD ^ ’ 


( 6 ) 


where D denotes the running UV cutoff, initially set by the width of the host band. Eq. ® 
yields a critical fixed point (SCR) at j* = r + O(r^), v* = 0, which separates the flows towards 
weak and strong coupling. Controlled calculations near SCR are therefore possible in a double 
expansion in r and j. Potential scattering is irrelevant at SCR and consequently does not play 
a role for leading critical exponents. 

Comparing these properties with the numerically determined phase diagrams one immedi¬ 
ately concludes that this critical theory only applies to 0 < r < rmaxj as otherwise no transition 
occurs at particle-hole symmetry. Indeed, a perturbative calculation of static critical proper¬ 
ties of SCR using the Kondo expansion shows excellent agreement with NRG results for small 
r < 0.2 HHj. 


3.4.2 SCR: Symmetric Anderson model 


The Anderson impurity model turns out to provide the relevant degrees of freedom to describe 
pseudogap Kondo criticality for all r > 0 [28]. To discuss the critical behaviour near SCR, 
we consider a symmetric Anderson model, Eq. ([^ with = — C//2, a momentun-independent 
hybridization v, and a particle-hole symmetric bath DOS as in Eq. ([^. The point Cd = U = 
n = 0 is referred to as the free-impurity fixed point (FImp), whereas the parameter sets v = 0 
and Ed = — C//2 = ±oo correspond to doubly degenerate local-moment states in the charge and 
spin channel, respectively. Therefore n = 0, = —oo can be identified with the LM fixed point, 

while u = 0, Cd = oo is dubbed LM’. 

Notably, the Anderson model is exactly solvable for any v a,t U = 0, known as resonant-level 
model. In the particle-hole symmetric case, its low-energy physics can be identified with that 
of the SSC fixed point introduced above: its properties correspond to a partial screening of the 
impurity degrees of freedom, with a residual entropy of S'imp = 2r In 2 [23EH]. 

A perturbative expansion is now possible in U around the SSC fixed point. The scaling 
dimension of the renormalized Coulomb interaction u at SSC is found to be dim)^] = —f = 
— (1 — 2r). The RG flow of u to two-loop order reads [33] 


du 

dluD 


= (1 — 2r) u 


3 (tt — 2 In 4) 3 


(7) 


This flow, together with the trivial flow near LM, LM’, and FImp is illustrated in Fig. [^ For 
all r > 0, LM is a stable fixed point, while SSG is stable only for r < rmax, as can be seen from 
Eq. Q. Therefore, a critical fixed point (SGR) emerges for 0 < r < rmax, Fig. [^, consistent 
with Fig. Its properties can now be accessed in a double expansion in f and u, and Eq. 0 
yields for the fixed-point coupling at = apjrryryiy?'- A perturbative calculation of static 
critical properties again yields excellent agreement with NRG results, here for r < Cmax [Ml- 
Owing to particle-hole symmetry, the behavior at e > 0 is formally identical to that at 
e < 0, Fig. § with the latter describing spin-Kondo physics while the former corresponds to 
charge-Kondo physics. SCR’ is a critical fixed point between a (partially) screened impurity at 
SSC and one with an unscreened charge degree of freedom at LM’. 


3.4.3 ACR: Asymmetric Anderson model. 

We now turn to the ACR fixed point present for r > r*. It was realized in Ref. m that 
the critical theory for ACR is that of a level crossing of a many-body singlet and a many- 
body doublet, minimally coupled to conduction electrons. Using the notation of Ref. [27128] . 













Will be inserted by the editor 


7 



(b) 1/2 < r < 1 



(c) r > 1 



Fig. 2. RG flow diagrams for the particle-hole symmetric Anderson model [21], in the plane spanned 
by the level energy e = — w/2 and the hybridization The symbols for stable fixed points correspond 
to those for the phases in Fig.[^ with LM and LM’ corresponding to local moments formed in the spin 
or charge channel; the flow of u near SSC is in Eq. 0- (a) 0 < r < 1/2: The critical fixed point SCR 
(SCR’) divides the flow to LM (LM’) from that to SSC. (b) 1/2 < r < 1: SCR and SCR’ merge with 
SSC as r — >■ 1/2“, such that SSC is now unstable, (c) r > 1: SSC merges with FImp at e = v = 0 
as r —>■ 1“. For all r > 1/2, LM and LM’ are the only stable phases in the presence of particle-hole 
symmetry. Reproduced from Ref. |29| . 


its Hamiltonian can be written as 


n = + £^o|cr)((7| +go [|cr)(s|c<,( 0 ) +h.c.] ( 8 ) 

k,(7 

where \a) = \ t),| i) and |s) represent the three allowed impurity states, eo is the tuning 
parameter (“mass”) of the QPT, i.e. the (bare) energy difference between doublet and singlet 
states. The QPT occurs at some Sq = Sc, with screening present for eo > £c- Remarkably, this 
theory is identical to a maximally particle-hole asymmetric Anderson impurity model, Eq. 0 : 
where the doubly occupied state has been projected out, U —>■ oo, and (eo,( 7 o) in Eq. ([^ have 
been identified with (€d,v) in Eq. 0. 

In this model, the point eo = 50 = 0 is dubbed valence-fluctuation fixed point (VEl). As 
above, go = 0, eo = —00 corresponds to LM, while go = 0, £0 = 00 describes a fully screened 
and particle-hole asymmetric singlet state, to be identified with ASC. 

A perturbative expansion is now possible in go around VFl. Power counting yields the scaling 
dimension of the renormalized hybridization dim[(/] = f = The one-loop flow equations for 

g and the renormalized mass e read 


dg 

dlnD 

de 

dlnD 


3 3 

-rg + ^g 

-e-g'^ + 3g‘^£ , 


(9) 


results to two-loop order can be found in Ref. [55]. The RG flow is shown in Fig. [^ - this 
flow has strong similarity to that of the standard Landau-Ginzburg model. The fact that g is 
relevant for r < 1 and irrelevant for r > 1 allows us to identify r = 1 as an upper critical 
dimension of the pseudogap Kondo problem, akin to d = 4 in the Landau-Ginzburg theory. For 
r < 1 a non-trivial fixed point (ACR) emerges at g*^ = |f and e* = — |r, Fig. Sk, similar to 
the celebrated Wilson-Fisher fixed point. Gritical properties, evaluated in a double expansion 
in f and g, again agree well with NRG results |2H]- In contrast, for r > 1 in Fig. [^, we have 
“Gaussian” behaviour controlled by the VFl fixed point, which here corresponds to a simple 
level crossing with corrections captured by plain perturbation theory in go- In the case r = 1, 
relevant to charge-neutral graphene, this perturbation theory is logarithmically divergent at 
criticality and needs to be resummed, as is standard at the upper critical dimension. 
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Fig. 3. RG flow for the asymmetric Anderson model |27I28| in the e-g^ plane, obtained from Eq. ([^ 
(a) r < 1: The critical fixed point ACR separates the flow towards LM from that to ASC. (b) r > 1: 
ASC merges with VFl as r —>■ 1 ~, which describes a level crossing with perturbative corrections. The 
behaviour near r = 1 is similar to that of the Landau-Ginzburg model near d = 4, with VFl and ASC 
corresponding to the Gaussian and Wilson-Fisher fixed points, respectively. Reproduced from Ref. [29]. 

3.5 Critical quasiparticles in the pseudogap Kondo model 

After having formally analyzed the critical field theories, we may now interpret the results. 
We start by discussing the single-particle electronic spectra, i.e., the question of critical quasi¬ 
particles. For fermionic impurity models, the relevant quantity is the conduction-electron T 
matrix; in the Anderson model ([^ this is related to the d-electron Green’s function according 
to T = v^Gd- 

In the stable LM and ACR phases of the pseudogap Kondo model, the T-matrix spec¬ 
tral density inherits the behavior of the bath, ImT oc \ujY |23EH|. In contrast, at criticality 
and for r < 1 the T matrix develops singular behavior, whose form is protected by the dia¬ 
grammatic structure of the T matrix (equivalent to a Ward identity [3S]) for both SCR and 
ACR fixed points: ImT oc |w|“’' |28l3fij . For r = 1 logarithmic corrections occur, leading to 
ImT oc 1/1 win cup, whereas a free-particle (5(cu) peak is restored for r > 1 [27]. We note that 
ImT cx |cu|“’’ is also realized in the SSC phase existing for r < 1/2. 

The interpretation of the critical behavior requires to distinguish the SCR and ACR fixed 
points. For SCR, both critical theories suggest that the criticality is exclusively carried by the 
spin degree of freedom of the impurity: The residual entropy is below In 2, and the Anderson- 
model RG flow is towards small charge fluctuations [28] . 

This is entirely different for ACR. Here, the residual entropy can exceed In 2, and for r close 
to or larger than unity the RG flow is towards large charge flutuations m- In particular, the 
“Gaussian” fixed point [i.e., the level-crossing description (|^] shows that spin and charge fluc¬ 
tuations are strongly coupled and equally important at criticality P7157] . Indeed, the criticality 
is carried by both fermionic and bosonic degrees of freedom, jcr) and js). The diagrammatic 
calculation shows that the spectrum of the physical d electrons arises from a convolution of 
critical bosons and fermions, which also explains how the exponent locking, ImT oc Iwj””, is 
compatible with non-trivially varying exponents for other observables. 


3.6 Pseudogap Kondo physics: Applications 

3.6.1 Kondo effect in graphene 

Graphene is the ideal host for the pseudogap Kondo problem: the local low-energy DOS follows 
Eq. ([^ with r = 1 at charge neutrality, and second-neighbor hopping t' 0 breaks particle- 
hole symmetry already on the level of the band structure. However, graphene is often doped, 
formally ^ ^ 0, due to the presence of a substrate or gates. Then, the DOS at the Fermi level is 
finite, and consequently a magnetic impurity described by the Kondo model ([l|) will be screened 
in the low-temperature limit for any value of the Kondo coupling Jq [38I39I40| . 

















Will be inserted by the editor 


9 


On general grounds, one expects that the presence of the quantum phase transition at 
charge neutrality also influences the behavior at finite /i. This has been studied in some detail 
in Ref. [40| using a combination of analytical and numerical renormalization-group techniques, 
and we summarize the most interesting aspects. Near the fj, = 0 phase transition, there exists 
a quantum critical regime in a plane spanned by the Kondo coupling and the doping. In this 
regime, the Kondo temperature is expected to follow the scaling prediction Tk = K±\fx\ for 
/i ^ 0. While this applies for r < 1, there are violations of scaling at the upper critical dimension 
r = 1, such that eventually Tk = K-|/i| continues to hold for p, < 0, while for /r > 0 logarithmic 
corrections and Kondo logarithms conspire such that Tk oc |/r|^ where a; « 2.6 is a universal 
exponent. Ref. [ID] has also given concrete predictions for the doping and parameter dependence 
of Tk away from criticality. 

Unfortunately, a clear-cut realization of Kondo physics in graphene is still missing. Exist¬ 
ing experiments have employed Co adatoms placed on top of a graphene sheet [41] as well as 
vacancies which are believed to create a localized magnetic moment [42143] . In the former case, 
individual Co atoms were investigated using STM techniques which indeed offer the unique 
opportunity to directly observe the critical quasiparticle spectra discussed above. However, the 
spectral features detected in Ref. [IT] mainly reflected charging effects and vibrational exci¬ 
tations. For the latter case of vacancies, signatures in low-temperature transport experiments 
have been interpreted in terms of Kondo physics |3D] , but alternative explanations of the data 
have been put forward as well [33]. 


3.6.2 Magnetic impurities on the surface of topological insulators 


The physics of magnetic moments coupled to boundary states of topological insulators has 
been discussed in numerous publications over the last years. For 2D topological insulators, 
the boundaries feature ID chiral edge channels which are expected to show Luttinger-liquid 
behavior. The corresponding Kondo problem has been argued to display both one-channel and 
two-channel Kondo phases [45] , with a possible quantum phase transition between the two |46| . 

In contrast, the surface states of strong 3D topological insulators admit a low-energy descrip¬ 
tion in terms of a 2d Dirac equation. The resulting electronic properties are therefore similar 
to that of graphene, with a few important differences: (i) there is a single Dirac cone (or, more 
generally, an odd number) per surface, and (ii) the role of the pseudospin (or sublattice) in the 
graphene case is taken by the physical spin, such that TI surface states display spin-momentum 
locking, and there is no additional spin degeneracy. The physics of Kondo impurities in this 
setting has been analyzed theoretically in a number of papers recently [47148149] . The main 
conclusion is that, despite the non-trivial topological structure of the TI surface states, the cor¬ 
responding local Kondo problem for a spin S' = 1/2 impurity can be mapped onto the standard 
pseudogap Kondo model of Sec. O Consequently, the physics is similar to Kondo screening in 
graphene, and concrete predictions for tunneling spectra and quasiparticle interference signals 
have been made [49] . To date, experimental investigations of isolated magnetic moments on TI 
surfaces at low temperatures are lacking. 


3.6.3 Quadratic band touching in 3D 

An interesting situation arises for quadratic band touching point^in 3D, as have been discussed 
for certain semiconductors such as HgTe m and also in the context in pyrochlore iridates Ell- 
Near a quadratic band touching point, the electronic DOS obeys p{uj) oc \oj\" with r = 1/2. 
Hence, one can expect that magnetic impurities in such a host realize a r = 1/2 pseudogap 
Kondo model (if the host electrons are assumed to be non-interacting). This Kondo model by 
itself features a quantum phase transition with strong scaling behavior and critical quasiparti¬ 
cles (provided that particle-hole symmetry is broken). It will be highly interesting to discuss 
concrete settings for this physics to emerge. 


2 


Quadratic band touching in 2D is realized, e.g., in bilayer graphene. 
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3.6.4 Weyl and Dirac semimetals in 3D 

Recently, some attention has been devoted to condensed-matter phases which feature linearly 
dispersing electrons in three space dimensions. This applies to Weyl semimetals |52j which 
may emerge, e.g., in certain pyrochlore iridates at intermediate correlation strength |53j . and 
to Dirac semimetals for which a number of materials realizations have been proposed |54l55j . 
For Dirac electrons in 3D, the DOS follows p{iS) oc |a;|’' with r = 2 if the chemical potential 
coincides with the energy of the Dirac point. 

As with graphene, this implies the existence of a quantum phase transition for Kondo 
impurities, i.e., screening will only occur for sufficiently large coupling and large particle-hole 
asymmetry. The corresponding quantum phase transition is, however, now above the upper 
critical dimension, i.e., it is a level crossing with perturbative corrections m- In particular, the 
residual local moment will jump across the transition, and the spectral function will display a 
<5 peak with subleading corrections, i.e., critical quasiparticles will not occur. 


3.6.5 Kondo effect in spin liquids 

A somewhat exotic realization of the Kondo problem arises in the context of magnetic impurities 
embedded in quantum magnets. If the host magnet is in a quantum spin liquid state, the host 
excitations are typically spin-1/2 spinons coupled to a U(l) and Z2 gauge field. The spinons 
couple to the impurity spin with a Kondo-like (i.e. four-point) interaction. The physics depends 
on the nature of the spinons, and a few cases have been discussed in the literature. 

Linearly dispersing bosonic spinons yield a rich phase diagram, with a variety of possible 
T = 0 phases, including the possibility of full Kondo screening, and quantum phase transitions 
[56| . In contrast, fermionic spinons lead to physics similar to standard Kondo expectations: In 
the presence of a spinon Fermi surface, the impurity spin gets always screened at low T [57] . In 
the case of 2D Dirac spinons of an algebraic spin liquid, a quantum phase transition not unlike 
that of the pseudogap Kondo problem emerges [58159160] . However, it should be noted that 
the influence of gauge fields beyond perturbation theory has been neglected in the published 
treatments. 


4 Critical quasiparticle theory for metallic QPT 

The conventional theory of continuous quantum phase transitions is phrased in terms of spatial 
and temporal fluctuations of an order parameter leading to a Ginzburg-Landau-Wilson (GLW) 
(j)^ field theory of the bosonic order parameter field. For metals this description has been pio¬ 
neered by Hertz [^, with refinements by Millis [62] and related work by Moriya (6^; we will 
below refer to this as HM theory. The order parameter field acts in deff = d + z dimensions, 
with d the spatial dimension and z the dynamical exponent. In case that deff > 4 , which is 
the upper critical dimension of held theory, the theory is Gaussian, i.e. the interaction of 
huctuations scales to zero at low energy under renormalization-group (RG) transformations. 

How do these predictions compare with observation? Quantum phase transitions have been 
studied experimentally for the last several decades. In particular metallic heavy fermion com¬ 
pounds show a broad variety of QPTs, mostly associated with some kind of antiferromagnetic 
order or also with superconductivity |3|. It was found that for some of the most intensely 
studied compounds, like GeCu6-a;Au2, (two-dimensional antiferromagnetic (AFM) fluctuations, 
dcff = 4) and YbRh2Si2 (three-dimensional AFM fluctuations, des = 5) HM theory does not 
seem to apply. One possible reason is that in these compounds one observes an apparently 
diverging quasiparticle effective mass m*, as a function of temperature T and excitation en¬ 
ergy u! , indicating critical behavior of the quasiparticle system. In such cases it is clear that 
a theory of critical phenomena should account for both, critical behavior of the bosonic and 
fermionic excitations. In the derivation of the GLW field theory for metallic systems fermions 
are integrated out, assuming they are non-critical. Such weak-coupling-based RG theories of 
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antiferromagnetic fluctuations in metals have failed to predict critical properties of the kind 
observed in the above-mentioned compounds [3]. 

One interesting option, suggested by experiment and available theory, is that the observed 
behavior is emerging at strong coupling. In the following a semiphenomenological theory of 
critical quasiparticles interacting with AFM fluctuations |65I66I67] will be described, resulting 
in a self-consistent relation for the quasiparticle effective mass. The latter allows for two types 
of solutions, depending on the initial condition at high energy, first the usual weak-coupling 
solution of HM theory, and secondly a new solution at strong coupling, involving critical quasi¬ 
particles with power law divergent effective mass. The crucial element of this theory is a Ward 
identity, relating the critical quasiparticle properties to the spin-quasiparticle interaction ver¬ 
tex, Xq oc m*/m , at finite momentum transfer q « Q, close to the momentum corresponding 
to the ordering wave vector Q. 


4.1 Critical quasiparticles 


Fermionic quasiparticle excitations are defined as poles of the single-particle Green’s function 

G 


G(k,w) 


1 

w — Ck — T’(k, ui) 


( 10 ) 


As we will see below the self-energy A'(k, w) depends only weakly on momentum k and we will 
neglect this dependence, A'(k, w) « XJ{u}). Expanding at small uj one may define the quasipar¬ 
ticle Green’s function G‘^^(k, w) 


G(k,a;)«—^^ = ZG^^’(k,a;), (11) 

oj — Ak + 

where the quasiparticle weight factor is defined as 

Z=[l-aReA/9w]-i = ^^, ( 12 ) 

m* (UJ) 

Ak = Z[ei^ + E(0)] is the renormalized quasiparticle energy and F = Zl-mE(uj) is the quasipar¬ 
ticle decay rate. We extend the usual Landau quasiparticle picture into the non-Fermi liquid 
regime by allowing Z to be energy dependent, Z = Z(uj) . Quasiparticles are well-defined as 
long as Z{uj) > 0 and F < |a;| , which is the case for a class of power law forms of the self-energy, 
E (X uj^~'^ , as long as 77 > 1/2 . To fully appreciate this statement one should observe that 
(1) the analyticity of ^( 0 ;) in the upper half plane leads to the relation ImA’(a;) oc Re 27(0;) if 
77 > 0 , (2) A = cr\uj\ for 0 < 77 < 1 , although Im27(a;) oc \uj\^~^, where cr < 1 if 77 > 1/2. . 


4.2 Critical spin and energy fluctuations 


The spectrum of antiferromagnetic spin fluctuations is given by m 


Imx(q,a;) 


NqujXq/j 

[r+eo{<i-Qn + {^>^Q/iy 


(13) 


where Nq is the quasiparticle density of states at the Fermi level, 7 « is a reference 
energy, and « 1/kp is the microscopic spin correlation length with vp, kp the Fermi velocity 
and wave vector. The control parameter r vanishes at the quantum critical point (QCP) as 
r (X (B — Bc)^" , where B is the field tuning the QCP (magnetic field, pressure, chemical 
composition) B^ marks the critical point and v is the correlation length exponent. The Landau 
damping term ujXq /y is renormalized by two factors of vertex functions Xq oc m*/m , appearing 
at both ends of a quasiparticle bubble diagram. 
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4.2.1 Vertex at large momentum 


The vertex correction Aq follows from the Ward identity connected with spin conservation 


d^p'duj' 
{2tt) 


d+l 


[f? - Aep>{Cl)]A{p,p'; Q) = n - Aep{Q) - Sip + Q/2) + Sip - Q/2) (14) 


where Aip,p'] Q) is the spin vertex function, Z\ep(Q) = ep+Q/ 2 ~ep-Q /2 • A singular behavior of 
the self-energy Sip) oc , 77 < 1 , implies Aq = limr 2 _>o /p, Aip,p'\Q) oc oc m*(a;)/m. 

The spectral weight given by (13) is peaked in momentum space at q = Q, implying that 


quasiparticles scattering off AFM spin fluctuations change their momentum by approximately 
Q , which means that if the initial quasiparticle state was on the Fermi surface, the final state 
will be “off-shell”, far from the Fermi surface, except for a small manifold of states on the Fermi 
surface, at so-called “hot spots”. 

There are at least two ways by which the critical scattering of quasiparticles may be dis¬ 
tributed all over the Fermi surface. The first one is by means of impurity scattering, which 
crudely speaking eliminates the constraints imposed by momentum conservation. This avenue 
has been explored in Refs. , where power-law critical behavior of the thermodynamic and 

transport properties has been derived, on the basis of a self-consistent equation for the quasipar¬ 
ticle effective mass. While the results for the critical exponents of that theory appeared to be in 
excellent agreement with experiment, the calculated prefactors of power laws were necessarily 
found to depend on impurity concentration, in contradiction with experiment. 


4.2.2 Critical behavior of tuning parameter 

The critical behavior of the tuning parameter r on the tuning field B (magnetic held, pressure, 
chemical composition), expressed in terms of the bare tuning parameter rg = (B — Be)/Be , 
is generated by the critical behavior of the dressed bubble diagram 7T(q, a;;B) in the static 
limit w = 0. The critical dependence of 7T on B can be found by considering the response 
of n to an inhnitesimal change of the tuning held, dB/dB, at hnite held B. The coupling 
vertex to the held is renormalized by a factor 1/Z, such that dB/dB oc 1/Z. Using the scaling 
behavior of Z with energy w as will be calculated below, Ziuj) oc and the general dynamical 
scaling relation uj oc where ^ is the spin correlation length, as well as the dehnition of 
the correlation length exponent ^ cx , one hnds dB/dB cx . Integrating with respect 
to B the scale-dependent contribution to B and hence to the control parameter r is found as 
r oc . On the other hand one has oc r oc r//' . Equating the two expressions for r one 

hnds 1 — rjzv = 2v from which follows v = 1/(2 -|- rjz). 


4.2.3 Energy fluctuations 


A second way by which quasiparticles may become critical over the whole Fermi surface is by 
scattering off-critical energy huctuations. The spin exchange energy operator is proportional 
to the scalar product of two nearest neighbor spin operators. The simultaneous propagation of 
two spin propagators with opposite momenta may therefore be viewed as energy propagation. 
The corresponding spectral function is given by m 

ImXE(q,w) = ^ G'fc+5iGfc+,i_Klmx(qi,a;i)Imx(qi - q,wi -a;)[6(a;i - w) - 6(a;i)], 
qi,^i 


where 6(a;) is the Bose function. The Green’s functions Gk+q^Gk+q^-q are off-shell for most 
values of the momentum qi and each may be replaced by l/ep- Performing the momentum 
integration by Fourier transform, we get 


ImX£;(q,w) 




1 


(15) 
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The possible relevance of these energy fluctuations near a QCP has been noted first by Hartnoll 
et al. |68j , in the context of a calculation of the dynamical conductivity in the perturbative RG 
regime. 


4.2.4 Dynamical scaling 


The spectra of spin or energy fluctuations at the critical point r = 0, as given by (131 and (151, 
have scaling properties with regard to scaling of energy and wavevector 


oj (X , z = 2/(1 — 2r]) . 


(16) 


It follows that the (bosonic) correlation length in the critical regime diverges as a function of 
temperature as ^ oc . 

Defining a correlation length 


we find that 


?(w;r) = [wAg/y + r] 


Imx(q,a;) oc 


ImX£;(q,a;) oc 


__ 

[r"(cc;r) + (q-Q)2p’ 

[r^(w;r) + g2](d+i)/2 ■ 


(17) 


(18) 


4.3 Critical effective mass 

The electron self-energy 27 may now be calculated in leading (one-loop) order. The imaginary 
part of 27 is obtained as 


/ I 

— ImG(k -b q, ct) 


'^)ImXi5(q, y)[b{y) + /(w - v)]. 


The interaction vertex \e = AgA«, where A„ is oc m*jm, as it arises through a Ward iden¬ 
tity connected to energy conservation, and Ag oc to*/to, as discussed above. The Fermi and 
Bose functions f{uj),b{uj) confine the ^/-integration at low T to the interval [0,a;]. The mo¬ 
mentum integration involves integration over the angle 9 enclosed by the momenta k, q. Using 
ImG(k -b q, w) « TrZd{uj — Zej,^q) and Ck-i-q ~ ek + Vk • q , where Vk • q =VFq cos 0 one finds 


J dcos9lmG(k,u}) 



vpq 



OJ -\- V 

Z 



(19) 


Here 0{x) is the step f unct ion. Performing the momentum and frequency integrations, using 
Imx£;(q, v) as given by (15) one finds for quasiparticles on the Fermi surface (ck = 0) 


Im i:(k,w) « 7AQ(wAy7)‘'+^/2C^(a;;r) 


( 20 ) 


where ^(w;r) has been defined by ( [T7| , observing that the step function in (15) is preempted 
by the condition q> , which defines the dominant regime of the q integral. 

In the critical regime (r = 0) we may now substitute the power law Ag(w) cx TO*(a;)/TO cx 
into (20) to obtain the scale-dependent contribution to Im E as 


ImS(k,w)(x 


(21) 
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The scale-dependent contribution to ReT’(w) follows from analyticity as ReT’(w) oc (AQ)^‘^+^(a;/ 7 )‘^ oc 
Substituting this into the definition of the effective mass one finds 
the self-consistency relation for m*{uj) m 


m*{uj) 

m 


1-b 




2d+l / X d-3/2 


( 22 ) 


This equation, for dimensions d > 3/2, has two types of solution, depending on the initial 
condition at high energy, at the upper end of the scaling regime. The first one is the usual 
H-M-model solution, for which the second term on the r.h.s. is always less than unity, implying 
that m*{uj)/m « 1 and the quasiparticle mass is not critical. If, however, the effective mass is 
already enhanced by some additional fluctuation contribution, such that the scale dependent 
term oc is already larger than unity, a new and self-consistent power law solution 

emerges. 


\ / \ 1/2—3/4d 

m V7/ 

_ 1 3 _ ri/4 d= 3 

2~ Ad~ \l /8 d=2 ' 


The critical exponents follow as 


2 _Ad 
^ ~ 1 - 2ry “ y ’ 

_ 1 _ 3 

2 “h Tfz 3 -|- 2(i 


(23) 


(24) 


4.3.1 Renormalization-group equation 


The self-consistent equation (22) for the effective mass describes scaling of the coupling constant 


Y = m*/m — 1 as a function of the scaling variable A = — ln(u;/ 7 ). The solution of (22) may 
be thought of as a result of solving a renormalization group (RG) equation for YilT). The 
corresponding RG equation is easily found as 


dA ’ Y-Yu 


(25) 


where Yu = l/2(i. The /3-function is seen to be negative for Y < Yu, so that for any initial 
condition at Aq, Y{Aq) < Yu the RG flow is directed away from Yu, towards the first stable 
fixed point at T = 0. Beyond the unstable fixed point at Yu the /3-function is positive, and 
the RG flow is directed away from Yu to the second stable fixed point at F —?> 00 . Since 
limy_>oo/3(F) = rfY, the power-law exponent describi ng the divergence of F is given by 7 , 
F cx , in agreement with the direct solution of Eq. (22). 


4.4 Thermodynamic scaling 


The free energy of the system is dominated by the critical contribution of the quasiparticles 
and may be calculated employing the expression for the entropy density S/V in terms of the 
self-energy 


= N(0) 


dw w(w — Re ^ 7 ( 0 ;)) 


27r 


T 2 cosh" ^ 


(26) 


For a critical discussion of (26) see [69]. Substituting the self-energy as determined by (20) 
using again the analytical properties of F'(w) in the complex frequency plane one finds S/V as 
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a function of the scaling variables temperature and control parameter Integrating 

with respect to temperature one finally gets the free energy density 

( 27 ) 

Here is the correlation length and Zf,Vf are the critical exponents of the fermionic degrees of 
freedom. One observes that the free energy density obeys hyperscaling. The spatial dimension of 
the fermions is d/ = 1, since momentum enters through the quasiparticle energy, which depends 
only on the distance from the Fermi surface. The dynamical exponent then follows from the 
quasiparticle Green’s function as Zf = 1/(1 — 77) , which determines oc. . Alternatively 

one may represent the free energy density in terms of the bosonic correlation length ^ cx 
in which case the correlation volume may be expressed as . We conclude that 

while the primary critical entities, the fermions, satisfy hyperscaling, the critical bosons do not. 
The bosons are driven by the fermions into a critical state intermediate between Gaussian and 
fully critical. We recall that in the cases considered here the bosonic fluctuations in the case 
of weak coupling to the fermions would be Gaussian. In the case of strong coupling we expect 
that the direct boson-boson interaction still flows to zero in the RG process. 

Employing the scaling form of the free energy, Eq. ( |27| , the critical contributions to the 
thermodynamic quantities, both in the critical regime (G) and in quantum disordered Eermi- 
liquid (EL) regime may be obtained. The specific heat coefficient 7 = C/T is found to diverge as 
7 oc T“’’ and 7 oc Tq 3 )/( 2 d+ 3 ) ^ regimes, respectively. The critical part of the 

polarization induced by the control parameter B (magnetization M, density p, ...) is found as 
M{T,Bc) — M{0,Bc) oc —T (C) and M(0,i?) —M(0, i?c) oc (FL). The corresponding 

susceptibilities x = dM/dB receive critical contributions x(r, He)—x(0i Be) oc _ 2 ’( 2 <i- 3 )/ 4 d 
and x(0, H)—x(0, Be) oc (FL). A quantity of special interest is the Griineisen ratio 

Bg = -{dMjdT)IC for which we find Bq oc C"! oc r-( 2 <i+ 3 )/ 4 d (c) and Bq = -Gr/{B - Be) 
(FL), where Gr = —{2d — 3)/(2d -b 3) is a universal coefficient [7^ . 


4.5 Scaling of transport properties 

The transport properties may be accessed in the framework of the critical quasiparticle picture 
by extending the relations known from Fermi liquid theory into the critical regime. The electrical 
resistivity p, for example, in FL theory is given by p = m*B/e^n, where e and n are the charge 


and the density of the carriers and B is the scattering relaxation rate. In the case that a 
typical momentum transfer in a scattering process is large, either intrinsically or by mediation 
of impurity scattering, one may approximate B by the imaginary part of the self-energy. Quite 
generally, one expects T, as a quantity of dimension energy, to satisfy the scaling relation 

B= r^$r(rr,roe'/")- 

(28) 

Gombining this with the scaling of the effective mass, one arrives at a 
inelastic part of the resistivity 

scaling form for the 

p{T,B)-p{T,Be) cx —r^<i>r{TC,roe^n 
m 

(29) 

1 ' j^{2d+3)/Ad ^(7^ 

Irg (FL) ' 

(30) 


4.6 Critical quasiparticles at quantum critical points in bulk metals: Summary 

The results presented in this section have been compared with experimental data for two heavy 
fermion compounds, GeCu 6 -a;Au 2 , [7S] and YbRh 2 Si 2 [7S]. In both cases the specific heat mea¬ 
surements show an apparentfy diverging effective mass in the limit T —0. In the first case the 
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critical point is located at concentration a; = 0.1 for ambient pressure and in zero magnetic 
field, or may be tuned by pressure and/or magnetic field away from x = 0.1. In the second case 
the transition is field tuned at ambient pressure and the critical field is Be = 60 mT. Neutron 
scattering studies indicate two-dimensional antiferromagnetic fluctuations in CeCug-xAua, |72j . 
with ferromagnetic fluctuations superposed m- Three-dimensional antiferromagnetic fluctua¬ 
tions have been identified in YbRh 2 Si 2 at low temperatures, T < 0.3 K, whereas ferromag¬ 
netic fluctuations prevail at higher T m- Applying the theory presented above assuming 
two-dimensional fluctuations (in the three-dimensional metal) the specific heat coefficient is 
found to diverge as 7 oc (rather than logarithmically, as conventionally assumed), in 

very good agreement with the data on CeCug-xAua, M- The theoretical result for the re¬ 
sistivity is p{T) — /9(0) cx (rather than linear), again in very good agreement with the 

data m- A crucial test of the theory is the scaling behavior of the dynamical structure fac¬ 
tor observed in inelastic neutron scattering studies m , which allowed an excellent fit to a 
power law scaling form with exponent of value a « 3/4, which is matched by our theory where 
a = l — 2ri = 3/4. Applying the theory assuming three-dimensional fluctuations to YbRh 2 Si 2 at 
low temperature and magnetic field, T < 0.3 K, R < 0.3 T one again finds excellent agreement 
with the available data [SlSSl- The specific heat coefficient is found to diverge as 7 oc T 
and 7 oc {B — in the C- and FL- regimes, respectively. The resistivity shows a sublinear 

T-dependence, p(T) — p(0) oc in the C-regime as seen in the data. The agreement of the 
scaling contributions with data on the magnetization, the spin susceptibility, the thermopower, 
and other quantities is also excellent [65166) . Of particular interest is the Griineisen ratio in the 
FL-regime, which is characterized by the universal number Gr , for which we find Gr = —1/3, 
while the experimental value is Gr ~ —0.3 . 

It appears that the self-consistent theory of quantum criticality [65166167) captures essential 
features of the quantum critical behavior of metals in situations where the quasiparticle effec¬ 
tive mass diverges, but the interaction of bosonic fluctuations is still weak. One is then in an 
intermediate regime of strong-coupling character with respect to the fermion-boson-interaction, 
but still in the weak-coupling regime of the boson-boson interaction. Anticipating that the one- 
loop approximation of the self-energy (employed in the above) is accurate, the theory would 
then be considered as a dynamically renormalized Gaussian theory. 


5 Conclusions 

We have reviewed the theoretical description of two cases of quantum phase transitions where 
fermionic spectral functions develop power-law singularities, a situation dubbed “critical quasi¬ 
particles” . The first case was about boundary phase transitions in the pseudogap Kondo model: 
Here the critical behavior is fully understood thanks to the existence of suitable epsilon expan¬ 
sion schemes. Possible experimental realizations include the Kondo effect in graphene, on the 
surface of topological insulators, and in systems with quadratic band-touching points. The sec¬ 
ond case was about an antiferromagnetic transition in a metal with strong coupling between 
fermions and spin fluctuations. Here the existing self-consistent theory is semiphenomenological, 
and a fully controlled derivation from a microscopic model is not yet available. 

We believe that various other quantum phase transitions, in particular those which involve 
the onset of Mott localization, should also display critical quasiparticles. Developing effective 
theories for those is an important task for future work. 
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by the Deutsche Forschungsgemeinschaft through the Research Unit FOR 960. 
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